Replica structure of one-dimensional disordered Ising models 
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We analyse the eigenvalue structure of the replicated transfer matrix of one-dimensional dis- 
ordered Ising models. In the limit of n — » replicas, an infinite sequence of transfer matrices is 
found, each corresponding to a different irreducible representation (labelled by a positive integer 
p) of the permutation group. We show that the free energy can be calculated from the replica 
symmetric subspace (p — 0). The other "replica symmetry broken" representations (p 7^ 0) are 
physically meaningful since their largest eigenvalues A^ control the disorder-averaged moments 
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In spite of extensive studies during the last twenty years, randomly disordered systems are not yet fully understood 



Pj . One of the main obstacles to the theoretical analysis of such systems is of course the lack of translational invariance 
of the Hamiltonian. The replica approach, which consists in replacing the original system and its randomly disordered 
Hamiltonian by the study of n (— > 0) identical and coupled systems with translationally invariant interactions, has 
been proven to be successful in providing a highly interesting mean-field theory of spin glasses (|||. The latter 
\Q ' relies on a fascinating, but mathematically obscur aspect of the replica mean-field theory, i.e. that the spin glass 
ON . transition coincides with the so-called replica symmetry breaking (RSB), or in other words, with the breaking of the 
permutation group symmetry of n — > elements. In this context, the crucial question is to find the correct scheme 
of breaking, allowing an analytical continuation of the theory when the number n of replicas tends to zero. To what 
extent the physical picture originated in the resolution of long-range models jlj applies to finite-dimensional systems 
and the occurrence itself of replica symmetry breaking are still under question |g| . 

In this note, we present some preliminary remarks about this important issue by analysing one-dimensional spin 
models with quenched random couplings and/or fields. As in the case of one-dimensional ferromagnetism, no transition 
is expected at finite temperature [|6|. However, since one-dimensional disordered systems may be solvable to a large 
extent with rigorous techniques, they constitute interesting examples on which the replica symmetry breaking approach 
can be tested before being applied to more realistic systems. 
■ To start with, let us consider the Hamiltonian 

C3 ' N 

H=-J2(JiS i S i+ i + h i S i ) (1) 
»=i 

with Ising spins Si = ±1 and periodic boundary conditions. The coupling constants and the external fields are 
randomly drawn with the probability distribution YiiLi "P{Ji)p{hi)- In general the 2 x 2-transfer matrices 

Ti(Ji,hi)= I e _ /3Ji _ j3ft . ^ph-phi J ( 2 ) 

at temperature 1/(3 are not simultaneously diagonalizable. The calculation of the thermodynamic properties, which 
requires the knowledge of the asymptotic properties of the product ]X Ti(Jj, hi), is rather involved but has been 
achieved for some particular choices of disorder distribution, using Dyson's method |||]-jl0| ■ 
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Hereafter, we shall follow a different route by noticing that, since the disorder is independently distributed from 
site to site, the free energy may be computed through the knowledge of the replicated and disorder-averaged transfer 
matrix T n =< T X (J, h)® n >, where < (.) > denotes / dJdh ■ V(J)p(h)[.). The entries of this 2™ x 2™-matrix are 



determined by the replicated spins {S a — ±1, a = 1, . . . , n}. We introduce the 2™ vectors \a\, et2, 



^ aj \/i ^ j, with up spins S a = +1 at and only at the sites a £ {01, 



They constitute an orthonormal 



basis of the underlying replicated space V n . The transfer matrix elements are now given by 



, a p \T n \bi, . 



/ n n 

, b a ) =< exp (3 J R a S a + (3hJ2R a 

\ r, — 1 1 




(3) 



with {R a }, {S b } corresponding to \ai, . . . ,a p ), \b\, . . . , b a ). They are obviously invariant under replica renum- 
bering, i.e. under all transformations of the permutation group S n given by its 2™-dimensional representation 
D(w)\ai, . . . ,a p ) = |7r(ai), . . . , 7r(a p )}, W G S n . Every irreducible decomposition of this representation D is iso- 
morphous to the most general eigenspace decomposition of V n . D already appears in the theory of atomic spectra, 
for a detailed presentation see p"l) . The number of up spins in the basis vectors of V n remains clearly invariant 
under replica permutations. So we find n + 1 subrepresentations A p which are carried by the spaces spanned by 



{I 



, a p ), 1 < ai < ... < a p < n}, p = 0, . 



But for p ^ 0, n these A p are still reducible. Consider e.g. 



the vector J2 a ^at 1°' ai ' ■ ■ ■ > a p-i); under permutations it behaves like \ai,. 



-1), so wc find A p = A p _i © D p 



Vp < n/2. Using in addition the symmetry p 

A 



n - p, 



± 



^F}, we obtain the complete decomposition of D: 
Ao = D G , Ai = D © Di, . . ., A p = D © D\ © . . . © -D m in(p,n-p), • • -i A„ = D , whose irreducibility has been proven 
in |nj. By changing the basis of V n with respect to this decomposition and taking at first the vectors of all Do-spaces, 
then those of all Di-spaces and so on, we can block-diagonalize the transfer matrix T n . 

Replica symmetry (RS) corresponds to the restriction to the first block since all Z?o~vectors are invariant under 
permutations. Due to the representation structure, the RS-transfer matrix has n + 1 non-degenerate eigenvalues and 
reads, cf. f| 



<t\ / n — a 



< exp{/3J(n + 4/x - 2r - 2a) + f3h(2a - n)} > 



(4) 



where /i_ = max(0, a + r — n) and p + = min(cr, r) and the indices a, r run from to n. The RS site-dependent 
partition function is given by the iterative description Z i+1 (r) = ^2& =0 Tn (cr,T)Zi(cr). Introducing the generating 
function Zi[x] = Zi(<r)x cr , the latter reads 



/•OO 

Z i+1 [x] = dy « e - phn {eP J + xe^ J ) n 8 (y - f(x)) » Z^y] 
Jo 



where 



-0j 
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(5) 
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In the thermodynamic limit, we call the right eigenfunction of T^^ which has the (maximal) eigenvalue unity 
and an integral normalized to one, 

(7) 



(8) 



$(x) = / dy < 5(x - f(y)) > $(y) . 
Jo 

The free energy density is given by the 0(n)-corrections in (^|), 

1 r°° 

/ = </!>--/ dx < log(e /3J + xe"' 3 ' 7 ) > $(x) 
P Jo 



Equations (0) and (||) for the random field Ising model are precisely the results derived inlijj by analysing the 
Lyapunov exponent of the infinite product of disordered transfer matrices Ti(Ji, hi) given in (|2). RS gives therefore 
the correct result for the free energy, and (|t]) is usually interpreted as the Dyson-Schmidt equation for the invariant 
density $ of the local fields || . A similar result was already established by Lin who showed the equivalence of 
an early replica method developped by Kac with Dyson's approach for the phonon spectrum of a chain of random 
masses and springs. 
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To be sure that replica symmetry is not violated, we have to check that the eigenvalue unity is not degenerate, 
i.e. reached in another eigenspace of T n ^Q. In the following, we shall therefore analyse the transfer matrix blocks 
corresponding to the non-trivial representations D p , p > 1. Each has n + 1 — 2p different eigenvalues of degeneracy 
(p) — (p-i) - By ordering the basis vectors according to their permutation properties one can achieve a further block- 
diagonalization of these blocks into (™) — ( "J identical blocks of size (n+1 — 2p) x (n + 1 — 2p) each one containing 
every eigenvalue of the D p -block exactly once. The transfer matrix blocks read 

TW (a, t) =« (2 sinh 2/3 J)" T^l 2p (a -p, T -p;J, h) » , (9) 

P < c,t < n — p, where T^(a, t; J, /i) is given by (|]) without averaging over the quenched disorder. In the limit 
n — > we obtain for every positive number p a different eigenvalue equation 

/>OC 

A (p) $(p) (x) = / dy « 0(3 _ /(y)) (/ ' (y)) p » $(p) {y) (10) 



where we again used the function f(x) defined in (||). For every eigenfunction &^(x) of the function -^^^(x) 

is an eigenfunction of the RS transfer matrix with the same eigenvalue. Only the largest eigenvalue (equal to one) 
of T^ , corresponding to the density (0), cannot be reached by this procedure. Therefore the largest eigenvalue of 
equals the second of and so on. 

To unveil the physical meaning of the maximal eigenvalues of all irreducible representations p > 1, one has to 
consider the disorder-averaged moments of the spin correlation functions \if =< ^ ({SiSj) — (Si)(Sj)) p More 
precisely, we shall see that they are dominated by 



(p) 
Xij « 



(*<")""* , (11) 



up to a multiplicative factor which depends on p. To show this, we consider the operator £( p ' = (E ® 1 — 1 ® £)® p 
where £ is the usual spin operator : Yi\S) — S\S) . The operator is defined on a set of at least 2p identical real 
replicas of the original disordered system. One can realize that it maps the spaces carrying Dq onto those carrying 
Dp. Since the ground state was found to be replica symmetric for any kind of disorder and is a priori non orthogonal 
to E^ (D p ), the correlations of two E^ p ) situated at sites i and j of the replicated Ising chain decay exponentially as 
(\(p))\j-i\. These correlations are in addition equal to the left-hand side of ([n]) as follows from the identity 

x\ P ? = 7T um Trace 

Ay 2P n^O 



4 P \T n y-% p) (T n ) N -i . (12) 



Therefore, the criterion \( p > — > 1 for replica symmetry breaking is equivalent to the divergence of the correlation length 
— (log A'-' 3 - 1 ) -1 of the p th moment of the connected two-point correlation function. The case p = 2 would correspond 
to the divergence of the spin glass susceptibility = 7^ j({SiSj) — (Si){Sj)) 2 , whereas the linear susceptibility 

X^ 1 ) = jf j({SiSj) — (Si )(Sj)) remains finite. This behaviour can be found in spin glass experiments as well as in 
mean- field theory Jjj. 

As a result of a replica calculation, the validity of (|l^,|ll|) is questionable. It is a simple exercise to check its 
correctness for the spin-glass chain without external fields : $( p )(a;) = S(x — 1) and A' p ) —<^. (tanh(/3J)) p 3>, for 
any p > 0, in agreement with known results jf|. For a more generic distribution of the disorder, exact calculations of 
the x^ p -*'s are much harder and, to our knowledge, have been performed for p = 1 only We show now that our 
replica calculation is exact in this case. Reproducing the approach of we compute the Fourier transform of x 
by adding a small sinusoidal field of wave- number q in the Hamiltonian( l]) and expanding the free-energy at second 
order in the field to obtain 

X (<Z) = dx da da du & q (x, a, a, u) <C 2u[ —r- — - (a 2 + a 2 ) —r- —tj 3> (13) 



K JJ + xe-P J J v ' \e? J + xe-0 J 

where the invariant measure $ g fulfills 

$>q(x, a, a,u) = J dy db db dv <& q (y,b : b : v) <C S(x — f{y)) 5(a — 2x — (bcosq + bsinq)f'(y)) 

1 



S(a - (bcosq- 6 sing)/' (j/)) S(u -2a + 2x- vf'(y) - ^(b 2 + b 2 )f"(y)) » . (14) 
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The susceptibility ( |l~3| ) depends on <I> g through the average values [tt] g (a;), [a 2 ] q (a;) and [a 2 ] q (x) only, where = 
J dadadu§ q (x, a, a, u)(.) (note that [l] 9 (a;) = 3>(x), see @)). These three functions obey three self-consistent equations 
obtained from whose integral kernels are regular. But their calculation still contains the averages [a] q (x) and 
[a] g (a;) fulfilling the complex equation 

POO 

/ dy (S(x -y)- « S(x - f(y)).f'(y) ») {[a] q (y) + i[a] q (y)) = 2x$(x) , (15) 
Jo 

cf. jl^). It becomes now clear that the pole Q of X^(o)i closest to the real axis, is reached when e l ® = sec (JTc|) . 
The asymptotic scalin g (Hll) coincides therefore with exactly known results for ■ It would be of high interest to 
extend the method of |14| (if possible) to test our replica predictions for larger values of p > 2. 

In this letter, we have studied the replica structure of one-dimensional Ising systems. From a technical point of view, 
we have first block-diagonalizcd the transfer matrix according to the irreducible representations of the permutation 
group of (finite) n elements and then sent n to zero in the eigenvalue equation for each representation. This procedure 
has allowed us to understand the physical meaning of the non trivial, "replica symmmetry broken" representations 
in terms of the moments of spin-spin correlations. Though one-dimensional Ising models do not display any replica 
symmetry breaking transition, we hope that the present approach will be of interest in future to understand if such 
a phase transition could occur in more realistic two or three-dimensional systems. 
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